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Abstract 

We present results from X-ray analysis of a Galactic middle-aged supernova rem- 
nant (SNR) G156.2+5.7 which is bright and largely extended in X-ray wavelengths, 
showing a clear circular shape (radius~50'). Using the Suzaku satellite, we observed 
this SNR in three pointings; partially covering the northwestern rim, the eastern rim, 
and the central portion of this SNR. In the northwestern rim and the central portion, 
we confirm that the X-ray spectra consist of soft and hard-tail emission, while in the 
eastern rim we find no significant hard-tail emission. The soft emission is well fitted 
by non-equilibrium ionization (NEI) model. In the central portion, a two-component 
(the interstellar medium and the metal-rich ejecta) NEI model fits the soft emission 
better than a one-component NEI model from a statistical point of view. The relative 
abundances in the ejecta component suggest that G156.2-I-5.7 is a remnant from a 
core-collapse SN explosion whose progenitor mass is less than 15 Mq. The origin of the 
hard-tail emission is highly likely non-thermal synchrotron emission from relativistic 
electrons. In the northwestern rim, the relativistic electrons seem to be accelerated 
by a forward shock with a slow velocity of ~500kmsec~^. 

Key words: ISM: abundances - ISM: individual (G156.2-I-5.7) - ISM: supernova 
remnants - X-rays: ISM 
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1. Introduction 



G156.2+5.7 is a large (radius~50'), X-ray bright supernova remnant (SNR) discovered by 
the ROSAT all-sky survey (Pfeffermann et al. 1991). The ROSAT 0. 1-2.4 keV band spectrum 
was well represented by emission from a thin thermal plasma with a temperature of about 
0.5 keV. A Sedov analysis based on the X-ray data showed a very low density medium of 
0.01 cm~^, an age of ~26000yr, and a distance to the remnant of ~3kpc. 

Subsequently, successful X-ray observations were performed with Ginga and ASCA. 
Using non- imaging instruments on board Ginga, Yamauchi et al. (1993) found that the scan 
profile in the 1.2-3 keV energy band was similar to the soft X-ray image observed with ROSAT, 
while the scan profile above 3keV revealed two bumps near to the north and south edges of 
the remnant. ASCA observed the northern and central portions of the remnant (Yamauchi et 
al. 1999). ASCA spectra from both north and central portions were characterized by soft and 
hard-tail components. The soft component was dominated by line emission, and originated 
from non-equilibrium ionization (NEl) plasma. A revised Sedov analysis based on the spectral 
parameters determined by the soft component resulted in a low ambient density of ~0.2cm~^, 
an age of ~15000yr, and a distance of ~1.3kpc. On the other hand, the hard-tail component 
was most likely represented by a power-law type spectrum, i.e., synchrotron emission from 
relativistic electrons. 

The distance to G156.2-I-5.7 has not yet been established. Gerardy & Fesen (2007) 
discovered considerable faint Ha line emission coincident with X-ray emission. Since some of 
them seems to be associated with interstellar clouds located at a distance of 0.3 kpc, they argued 
that the distance to the remnant might be closer than that was estimated previously. We here 
take Ikpc as a distance to the remnant, following the most recent paper on G156.2-I-5.7 (Xu 
et al. 2007). 

Radio observations show a bilateral shell structure having relatively weak emission at 
the central portion of the remnant (Reich et al. 1992). Although the radio surface brightness 
was the lowest among all known SNRs, a magnetic field in the shell was estimated to be 16 /zG 
(at a distance of Ikpc) which is just slightly higher than a typical interstellar magnetic field. 
Therefore, Reich et al. (1992) concluded that the low surface brightness was due to the low 
ambient density. Polarized radio emission shows a nearly perfect tangential magnetic field 
configuration along the shell, which is typical for evolved shell- type SNRs (Reich et al. 1992; 
Xu et al. 2007). All the observations agree that G156.2-I-5.7 is a middle-aged SNR whose SN 
explosion occurred in a low-density ambient medium. 

We focus on both thermal and non-thermal X-ray emission in this SNR, using the most 
recent X-ray astronomy satellite Suzaku (Mitsuda et al. 2007). The fact that non-thermal X-ray 
emission is detected from this middle-aged SNR is very interesting, since a currently accepted 
theory of diffusive shock acceleration requires high shock velocity of at least ~2000kmsec~^ 
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Fig. 1. ROSAT PSPC image of t lie entire G156.2+5.7. The energy range is 0.05-3. OkeV. The effects of 
vignetting and exposure are corrected. The image is binned by 15" and is smoothed by a Gaussian kernel 
of <T = 45". The intensity scale is square root. The Suzaku/XIS FOV (17'.8x 17'.8) are shown as white 
boxes. 

which is not expected in middle-aged SNRs to produce non-thermal emission in the X-ray 
domain. Fortunately, the thermal emission from this SNR allows us to investigate ambient 
medium densities as well as the shock velocities which are important to study the theory of 
diffusive shock acceleration. Another interesting aspect of this remnant is a clear evidence of 
Si and S rich ejecta detected in the central portion of the remnant with ASCA (Yamauchi et 
al. 1999). We can investigate the progenitor star which produced this SNR, by comparing the 
relative metal abundances of the ejecta with those expected in nucleosynthetic models. 

2. Observations and Data Screening 

We observed the eastern (E) rim, the northwestern (NW) rim, and the central portion 
of G156.2+5.7 during 2007 February 16-20 with the Suzaku satellite. We here concentrate 
on the data taken by the X-ray Imaging Spectrometer (XIS; Koyama et al. 2007) onboard 
Suzaku. The XIS consists of two front-illuminated (FI; XISO and XIS3) CCD cameras and one 
back-illuminated (BI; XISl) CCD camera. Each camera covers an identical imaging area of 
17'.8x 17'.8. The fields of view (FOV) of XIS for the three observations are shown as white boxes 
in the entire X-ray image of the SNR obtained by the ROSAT PSPC (figure 1). We reprocessed 
the XIS data, using the latest CTI calibration file of version 20080131. We cleaned the repro- 
cessed data by the standard criteria^ recommended by the calibration team of Suzaku/XIS. 

^ See the Suzaku Data Reduction Manual which can be found from 
http:/ /heasarc. gsfc.nasa.gov/docs/suzaku/analysis/abc. 
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Fig. 2. NXB-subtracted Suzaku XISO images of the three pointmgs. The energy range is 0.4-8.0 keV. 
The effects of exposure, vignetting, and contamination are corrected. The image is binned by 8" and 
is smoothed by a Gaussian kernel of ct = 24". The intensity scale is square root. White curves are the 
regions where we extract spectra. We exclude two point sources enclosed by dashed circles from our 
spectral analyses. 

After the screening, the remaining exposure times were 53.3 ks for the E rim, 50.5 ksec for the 
NW rim, and 51.2 ksec for the central portion, respectively. 

Figure 2 shows the vignetting-corrected XISO images of the three observations. Two 
corners of each FOV where the calibration source of ^^Fe is illuminated are masked in these 
images. 

3. Spectral Analysis and Results 

Figure 3 shows XISO spectra extracted from the entire FOV for each observation. The 
non X-ray background (NXB; see the next subsection for an explanation) spectra are also shown 
as red crosses. The signal-to-noise (black-to-red in figure 3) ratios in an energy band of 3-5 keV 



4 



^ 0.1 



o 0.01 
u 



10-: 



S V 



(a) East 



!-+++-H- 



0.5 



2 5 
Energy (keV) 



10 



=S 0.1 



o 0.01 

u 



10- 



(b) Northwest 



f-+4.-H-"' 



0.5 1 2 5 

Energy (keV) 



10 



^ 0.1 



S 0.01 

o 



10- 



, A, ] 



(c) Center 



'A ' s 



0.5 1 2 5 

Energy (keV) 



10 



Fig. 3. XISO spectra from the entire FOV (blaclc crosses), (a), (b), and (c) arc responsible for tlic E rim, 
the NW rim, and the central portion, respectively. Non X-ray background spectra are also shown as red 



significantly vary among the three FOV, indicating a non-uniformity of the hard-tail emission 
detected by Ginga and ASCA (Yamauchi et al. 1993; 1999). Also notable is a variation of 
equivalent widths of K-shell lines from Si (at ~1.85keV) and S (at ~2.4keV); those in the 
central spectrum are much larger than those in the rim spectra. This implies the existence of 
the fossil ejecta in the central portion as was detected with ASCA (Yamauchi et al. 1999). In 
the following spectral fitting procedure, we use photons in the energy range of 0.4-8.0 keV, and 
employ the XSPEC (version 12.4.0) software. 

3.1. Background 

As background, we consider three kinds of sources; NXB caused by charged particles 
and 7-rays hitting the detectors, cosmic X-ray background (CXB) which accounts for emission 
from unresolved point sources, and local hot bubble (LHB) that is a hot gas surrounding the 
earth. Note that Galactic ridge X-ray emission (GRXE) is negligible for this SNR; based 
on literature by Sugizaki et al. (2001) and Revnivtsev et al. (2006), we estimate the flux 
of GRXE at the Galactic longitude, /, / = 100° within Galactic latitudes, b, \b\ < 1° to be 
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4xl0~^^ergcm~^ sec~^ arcmin"^ in a range of 0.7-10 keV, which is about 0.3% of that from 
G156.2+5.7 (Yamauchi et al. 1999). Furthermore, the larger values of /(= 156.2) and b{= 5.7) 
for this SNR would prefer the smaller flux of GRXE. We generate NXB spectra suitable for 
our observations, by employing xisnxbgen provided by the XIS team (Tawa et al. 2008). The 
generated NXB spectra are shown in figure 3, in which we can see that the NXB spectra above 
8keV, where emission from astrophysical sources is negligible, match the data. Therefore, we 
are confident that NXB emission is properly subtracted by the NXB spectra generated by 
xisnxbgen. We subtract NXB emission from the same area as the source area in detector 
coordinates. For the CXB emission, we employ broken power-law model with photon indices F 
(<0.7keV) = 2.0 and F (>0.7keV) = 1.4 (e.g.. Miller et al. 2008). As for the LHB component, 
we use the APEC model (Smith et al. 2001) with solar composition (Anders & Grevesse 1989) 
and an electron temperature of 0.1 keV (e.g., Tawa 2008). We estimate the emission measure of 
the LHB component around G156.2+5.7 to be 2xl0~^pccm~^, based on the surface brightness 
in the ROSAT PSPC Rl and R2 bands of ~2.8 counts sec~^ arcmin"^. We exclude two point- 
like sources visible in our FOV from our spectral analyses; the regions excluded are enclosed 
by dashed circles as shown in Fig. 2. 

3.2. Emm 

We perform annular spectral analyses for the E and NW rims. We choose the center 
of the X-ray SNR at position (RA, DEC) = (4'^59'^05% 51°50'57") [J2000.0] determined by a 
ROSAT PSPC surface brightness contour. For the E rim, the annular regions are spaced by 
4' from a radial position, r = 48' to r = 32', resulting in four annular regions (from El at the 
outermost region to E4 at the innermost region). These annular regions are shown in figure 2 
as white curves. 

We perform spectral fittings for the NXB-subtracted spectra. Since previous ASCA 
observations showed non-equilibrium ionization (NEI) states in this SNR, we apply an absorbed 
NEI model (VNEI model [e.g., Borkowski et al. 2001]) having LHB and CXB components 
for all the four spectra. The free parameters are the hydrogen column density, A'h; electron 
temperature, kTg; the ionization timescale, r (r is the product of the electron density multiplied 
by the time after the shock heating); the emission measure, EM (EM= / nf^nndl, where Ue and 
TT-H are the number densities of electrons and protons, respectively and dl is the plasma depth); 
abundances of N, O, Ne, Mg, Si, S, and Fe. The other elemental abundances are assumed to 
the solar values (Anders & Grevesse 1989), since emission lines from these elements cannot be 
seen in the Suzaku spectra. We first allow the normalization of the CXB component to vary. 
We then find that the mean flux of the CXB component derived in the four annular regions 
is 4.7x10^^^ ergcm"^ sec~^ arcmin"^ in 2.0-10.0 keV. This is consistent with the CXB flux of 
4.4x10"^^ ergcm"^ sec~^ arcmin"^ estimated by blank sky observations with Suzaku (Tawa 
2008). Since the CXB flux fluctuates by about 10% from location to location, we take the 
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mean CXB flux derived in the four annular regions as a local flux around this SNR. After fixing 
the normalization of the CXB to the mean flux, we fit all the spectra with a one-component 
VNEI model. This model gives us fairly good fits for all the spectra; the reduced-x^ values 
(d.o.f.) range from 1.04 (753) to 1.27 (837). No additional power-law component for possible 
hard-tail emission is required at a 99% confidence level. Therefore, we consider that this model 
well represents the spectra in the E rim. We use this CXB flux value in the following analysis. 
Figure 4 (a) shows the example spectrum (from E3) with the best-fit model. The best-fit 
parameters for the four spectra are listed in Table 1. 

We find that the values of kT^ significantly decreases toward the shock front from 
0.44 keV to 0.28 keV. This trend is generally seen either in Sedov-phase SNRs or behind the 
shock fronts interacting with dense interstellar medium (ISM) (e.g., Miyata & Tsunemi 1999; 
Levenson et al. 2005; Hwang et al. 2005). The measured abundances are almost constant in 
the FOV. The typical abundance of 0.5 times the solar value is consistent with those derived 
in the northern rim of this SNR obtained with ASCA (Yamauchi et al. 1999). The sub solar 
abundance leads us to consider that the origin of the plasma in the E rim is the ISM swept-up 
by the forward shock. We should note that the N abundance is depleted to at least 0.1 times 
the solar value which is not expected in the ISM plasma. In our data, it is quite difficult to 
establish whether or not this strong depletion is real, given the fact that we do not know the 
composition of contaminants built up on the optical blocking filters of the XIS (Koyama et al. 
2007) in addition to weak N emission lines from G156.2+5.7. The value of A^h is almost con- 
stant in the FOV, and is consistent with those derived by the ASCA observations (Yamauchi 
et al. 1999). 

3.3. NW Rim 

In the NW rim, we also extract four spectra from annular regions (from NWl at the 
outermost region to NW4 at the innermost region). The width of the innermost region is 
6', while those of the other regions are 4'. These annular regions cover from r = 30' to r 
= 48', and are shown in figure 2 as white curves. We first apply the same model as we used 
for the spectral modeling in the E rim, i.e., a one-component VNEI model for the emission 
from G156.2+5.7. This model gives us an acceptable fit for the innermost region (NW4). 
Therefore, we consider this model as the best-representative model for the NW4 region. The 
best-fit parameters are listed in Table 1. On the other hand, the fit levels for the other three 
spectra are not acceptable; we see excess emission from the best-fit models above 3keV. This 
is at least qualitatively consistent with the fact that Yamauchi et al. (1999) detected hard-tail 
emission above 3keV in the ASCA spectra. Therefore, we add either a thermal (NEI) or a 
power-law component for the hard-tail emission. We set metal abundances of the additional 
NEI component to 0.5 times the solar value which is a typical abundance of the ISM around 
G156.2+5.7. The kT^, r, and EM for the additional NEI component, or photon index and 
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Fig. 4. Example spectra with their best-fit models. The best-fit curves are shown as solid lines (total 
components) or dotted lines (individual components). The lower panels show the residuals. The spectral 
extraction regions are (a) from E3 in figure 2 (a), (b) from NW2 in figure 2 (b), and (cl) or (c2) from the 
ellipse region or the surrounding region in figure 2 (c). 

normalization for the additional power-law component are treated as free parameters. Both 
VNEI-I-NEI and VNEI-|-power-law model dramatically improve the fits for all the three spectra; 
e.g., the value of x^/d.o.f. for the NW2 region is reduced from 1240/911 to 1065/909. The 
example spectrum with the best-fit VNEI-|-power-law model for the NWS region is shown in 
figure 4 (b). The best-fit parameters for the three spectra are shown in Table 2. 

The surface brightness of the hard-tail component, ~5xl0~^^ ergcm"^ sec"^ arcmin"^, 
is consistent with that detected with ASCA/GIS in the northern portion of the remnant 
(Yamauchi et al. 1999). We find that the VNEI-|-power-law model fits the data slightly bet- 
ter than the VNEI-I-NEI model. However, it is not very safe to exclude the possibility of the 
VNEI-I-NEI model (i.e., the thermal origin for the hard-tail emission) from a statistical point 
of view. Thus, we consider about astrophysical aspects. If thermal origin is the case, the high 
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Table 1. Spectral-fit parameters for the E and NW rim. 



Parameter 










TT A 
hj4 


IN W4 




3.7±0.1 


3.2±0.1 


2.5±0.1 


3.4±0.1 


4.1±0.1 


A:To[keV] 




0.28±0.(Jl 


A O C 1 A A 1 

0.35±0.U1 


A /I c: 1 A A1 

U.45±0.01 


A A A 1 A A 1 


A c: 1 1 A A 1 

U.51±0.01 


N 




All 1 A A C 

0.11±0.05 


<0.0d 


<0.02 


<0.03 


^ A f\ A 

<0.04 







A O O 1 A A 1 


A AO 1 A A1 

0.43±0.01 


A O O 1 A /^^ 1 

0.38±0.01 


A OO 1 A A1 

0.32±0.01 


A OO 1 A AO 


Ne 




A /^O 1 A AO 


A CC\ 1 A AO 


A fTi 1 A AO 

0.d0±0.02 


A /I A 1 A AO 

0.4U±0.02 


A /I O 1 A AO 

0.43±0.03 


Mg 




A P" T 1 A A P" 

0.57±0.05 


('\ r ^ 1 A A 1 

0.57±0.01 


A /t A 1 (~\ ('\ A 

0.49±0.04 


A O O 1 ('\ r\ A 

0.32±0.04 


A 1 1 A A /I 

0.41±0.04 


Si 




0.8±0.2 


0.5±0.1 


0.5±0.1 


0.4±0.1 


0.9±0.1 


S 




<1 


<0.5 


<1 


<0.7 


1.0±0.7 


Fe 




n /I c no 
U.45±U.Uz; 


U.51±L).01 


0.4o±0.02 


n o /I no 
0.34±0.0z 


A oo^n no 
O.o8±0.02 


log(r/cm~'^ 


sec) . . . 


10.96±0.03 


10.80±0.03 


10.62±0.03 


10.50±0.03 


10.53±0.04 


EMt[xlOi« 


cm^^] . 


3.26±0.07 


1.63±0.03 


0.72±0.02 


1.22±0.04 


0.74±0.02 


Surface bri^ 


jhtness^ 


3.0x10"" 


2.4x10-1^^ 


1.3x10"" 


1.9x10"" 


1.3x10"" 


xVd.o.f. .. 




781/753 


1065/837 


886/763 


684/569 


632/603 



* Other elements are fixed to those of solar values. 

The values of abundances are multiples of solar value. 

The errors are in the range Ax^ < 2.7 on one parameter. 

tEM denotes the emission measure ^ ricrijidl. 

^In unit of ergcm^'^ sec^^ arcmin^^ in a range of 0.5— lOkeV. 



kTf, of above 3keV suggests that G156.2+5.7 is as young as historical young SNRs. As aheady 
pointed out by Gerardy & Fesen (2007), this is unhkely, since there have been no historical 
record of the SN explosion toward G156.2+5.7. Furthermore, the spatial variation of /cTg, i.e., 
the value of kTc is at least two times higher at the outermost region than those in the inner 
regions (see. Table 2), conflicts what we expect behind shock fronts; kT^ generally decreases 
toward the shock front due to deceleration of the shock as well as the effect of thermal non- 
equilibration between ions and electrons. On the other hand, if non-thermal origin is the case, 
the spectral flattening towards the outer region (see. Table 2) is consistent with that seen at 
the rim of Tycho's SNR (Cassam-Cheinai et al. 2007). Therefore, we believe that the origin of 
the hard-tail emission is non-thermal emission rather than thermal emission. 

The abundances derived by the VNEI component responsible for the soft emission from 
G156.2+5.7 are similar to those observed in the E rim, suggesting that the origin of the plasma 
in the NW rim is also the swept-up ISM. The values of kT^ are also similar to those in the E 
rim. This fact implies similar forward shock velocities in both the NW rim and the E rim. 

3.4. Center 

We extract two spectra in the central FOV. One is a central ellipse region, where we can 
see enhanced surface brightness relative to the surrounding region in the ROSAT PSPC image 
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Table 2. Spcctral-fit parameters for the NW rim and the central portion. 



U 111 1.^ w lH_j IILj 


Tf! TIIOIT'I' 

1_ OjI (XllHj Lj'^I 


NWl 


NW2 


NW3 


In'.lli-ncri 


Sill I'TTimi n m 

kJ IJ-l 1 \J LlllVlllltl 






IVTnHpl ■ 


VNF,T-I-NFJ 










JVH['^i^ CIIl J ... 


q 7_|_n 1 


comu. 1 


o.ymu. i 


4 1 -1-0 1 


q Q+o 1 


VNTTT 


kT \h-n\n 


u.ozinu.ui 


n zio+0 01 


n 4Q-|-n 01 


n 4fi+0 01 


n c;9-i-o 01 




1\T 

iy 




^0 01 


.r'O 09 


^0 09 


^0 01 




\j 


n Q/l-l-n 09 


n qt^+o 01 


n q9-i-o 09 


n q9+o 01 


n 98-1-0 01 




1 N C 


n c;9+o 09 


n fic;+o 09 


n ^0-1-0 09 


n qo+o 09 


44-1-0 09 




A/lrr 

iVig 


4/1+0 0^^ 


n fiq+o 04 


n 4fi-l-0 0"? 


n q7+o 09 


n 4f;-l-0 0"? 




ol 


n f;-i-n 9 


7-1-0 1 
U. / ztU. 1 


R-l-0 1 


9 1 -1-0 1 
Z. iztU. i 


"? 0+0 1 




q 

o 


<^o s 


^0 4 


^1 9 




0.0+u.u 




Fp 


n qy+fl 09 


n 47+0 09 


Q qa+o 01 


n c.o+0 01 

W. ?JU— 1— U.VJ X 


Q c;q+o 01 




iOgt^T/ CIIl bcC J . . . . 


lU.oZZCU.U'i 


1 n 70-1-0 0"? 


1 n f^/L-l-fl OA 


lu.oomu.uz 


1 71 -1-0 09 
lU. ( IHZU.UZ 




pi\/rt [ V 1 ni8f,m-5i 

H/iVi 1 AlU CIIl 1 . . 


Z.UOHzU.UD 


1 1 4+n 09 


1 1 7-1-0 V\'X 

1.1/ ZlZU.UO 


1 91 09 
1 .ZlUZU.UZ 


fi'^-l-0 01 




Surface brightness"^ 


Z.O X iU 


1 7n/10^13 
1. ( X lU 


1. / X lU 


9 1 N/ 1 n~i3 
z . 1 X lU 


1 1 1 0~13 
1 . 1 X lU 








q 9+0.4 
"^■^-0.5 




e+4 
0-2 






Metal abundance . 


0.5 (fixed) 


0.5 (fixed) 


0.5 (fixed) 


0.5 (fixed) 


0.5 (fixed) 




log(T/cm~^ sec) .... 


>12 


>12 


>12 


>12 


>12 




EMt[xlOi^ cni^S] . 


0.030±0.005 


0.06±0.01 


0.053±0.006 


0.027±0.004 


0.014±0.003 




Surface brightness^ 


3.7x10-15 


5.2x10-15 


4.8x10-15 


3.1x10-15 


1.8x10-15 


xVd.o.f 




676/588 


1065/909 


981/796 


1375/1101 


1796/1332 






Model: VNEI+powcr-law 


Absorption 


iVH[xl02icin-2] ... 


3.9±0.1 


3.4±0.1 


3.9±0.1 


4.0±0.1 


3.9±0.1 


VNEI 


A:To[keV] 


0.31±0.01 


0.37±0.01 


0.43±0.01 


0.48±0.01 


0.50±0.01 




N 


<0.04 


<0.02 


<0.04 


<0.03 


<0.03 







0.36±0.02 


0.43±0.02 


0.40±0.02 


0.36±0.01 


0.38±0.02 




Ne 


0.53±0.02 


0.85±0.03 


0.66±0.03 


0.46±0.02 


0.64±0.03 




Mg 


0.44±0.05 


0.80±0.05 


0.59±0.05 


0.44±0.03 


0.64±0.04 




Si 


0.5±0.2 


0.8±0.2 


0.8±0.1 


2.4±0.1 


4.2±0.2 




S 


<1 


<0.7 


<2 


2.1±0.5 


5.9±0.8 




Fe 


0.37±0.02 


0.60±0.02 


0.49±0.02 


0.58±0.01 


0.80±0.02 




log(r/cm^'^sec) .... 


10.84±0.04 


10.79±0.03 


10.67±0.04 


10.58±0.02 


10.74±0.02 




EMt[xlOis cni-5] . 


2.25±0.06 


1.05±0.02 


0.90±0.02 


1.00±0.02 


0.49±0.01 




Surface brightness^ 


2.3x10-13 


1.8x10-13 


1.6x10-13 


2.0x10-13 


1.1x10-13 


Power-law 


Photon Index 


1 7+0.2 
^- ' -0.1 




2.6±0.1 


2.3±0.1 


2.6±0.2 




Norm-'- 


7±1 


38±3 


27±2 


13±2 


13±2 




Surface brightness^ 


3.9x10-15 


9.9x10-15 


7.9x10-15 


4.5x10-15 


3.5x10-15 


xVd.o.f 




674/589 


1036/909 


954/796 


1375/1101 


1781/1332 







* Other elements are fixed to those of solar values. 
The values of abundances are multiples of solar value. 
The errors are in the range Ax^ < 2.7 on one parameter. 
tEM denotes the emission measure f nrnwdl. 



(see, figure 1), the other is the rest of the FOV. We first apply the same one-component VNEI 
model as we did for the E rim spectra. Then, we obtain the reduced-x^ values (d.o.f.) for 
the central ellipse and the surrounding region to be 1.33 (1103) and 1.38 (1334), respectively. 
Taking into consideration that systematic uncertainties of the response, the fit level might be 
acceptable from statistical point of view. However, apparent excess emission in the hard energy 
band (above 3keV) is found. Therefore, we add either an NEI or a power-law component for 
the hard-tail emission. The best-fit parameters are shown in Table 2. Although the additional 
component (either NEI or power-law) significantly improves the fits for both spectra (reduced- 
values are 1.25 for the ellipse region and 1.34 for the surrounding region), we still see obvious 
discrepancy between our data and the best-fit models around 0.5-0.6 keV (He-like O Ka line) 
and 1.2-1.3 keV (Fe L). This means that the VNEI-|-(NEI or power-law) model is too simple to 
reproduce the soft (thermal) emission from G156.2-I-5.7. 

Since Yamauchi et al. (1999) found metal-rich ejecta at the central portion of the SNR, 
it is reasonable to consider that at least two kinds of plasmas (the swept-up ISM and the 
metal-rich ejecta) exist along the line of sight at the center of the SNR. We thus apply a two- 
component (ISM and ejecta component) VNEI model to represent the thermal emission from 
G156.2-I-5.7. Assuming that abundances of the ISM are uniform around G156.2-I-5.7, we set 
abundances of the ISM component to the typical values derived in the E or NW rims (N=0.05, 
0=0.4, Ne=0.6, Mg=0.5, Si=0.6, S=0.6, and Fe=0.5 x the solar values, and the other metals 
are assumed to be the solar values). In order to reduce the number of free parameters and 
obtain meaningful results, we also set the value of kTc for the ISM component to 0.3 keV that 
is typical for the outermost regions in the E and NW rim. The kT^ for the ejecta component 
is left as a free parameter. The r and EM are left as free parameters for both ISM and ejecta 
components. We also left the metal abundances of O, Ne, Mg, Si, S, and Fe for the ejecta 
component as free parameters. Since it is quite difficult to constrain the abundance of N due 
to poor statistics below 0.5 keV for the ejecta component, we set the abundance of N to be 
equal to that of O. In this way, we apply the 2VNEI-|-(NEI or power-law) model for the two 
spectra extracted in the central FOV. This model significantly improves the fits. The best-fit 
parameters are listed in Table 3. The example spectrum from the surrounding region is shown 
in figure 4 (c). Since this model gives us nearly acceptable fits for both spectra and there is no 
strong reason to add further components, we stop introducing additional components. 

If we assume that the hard-tail component is represented by the thermal model, the 
value of kTc is measured to be unreasonably high (above ~6keV). Therefore, we consider that 
the hard-tail emission is most likely explained by the power-law type spectrum rather than the 
thermal spectrum. The metal abundances for the ejecta component are indeed higher than 
the solar values, confirming that it originates from the metal-rich ejecta. The abundances in 
the ejecta component are significantly different between the central ellipse and the surrounding 
region; the abundances for the central ellipse are all lower than those for the surrounding region. 
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The relatively low abundances at the central ellipse might be caused by either a mixing of the 
ISM into the ejecta component or technical difficulty in separating the ISM component and the 
eject a component. 

4. Discussion 

We detected line-dominated soft emission in all the three pointings i.e., the E rim, the 
NW rim, and the central portion of this SNR. In the E and NW rims, the spectra of the soft 
emission show similar shapes, and are represented by thermal emission originating from the 
swept-up ISM. On the other hand, in the central portion, the spectra of the soft emission show 
an additional component corresponding to metal-rich ejecta. In addition to the soft emission, 
we detected hard-tail emission in the NW rim and the central portion. The origin of this hard- 
tail emission is more consistent with non-thermal synchrotron emission rather than thermal 
emission. 

4.I. Thermal Emission 

Our annular analyses for the rim regions revealed spatial variation of spectral parameters. 
The /cTg decreases from the innermost region (~0.5keV) toward the shock front (~0.3keV), 
while the EM increases toward the shock front. These trends are at least qualitatively consistent 
with those expected in the Sedov-phase SNRs. Then, we compare the EM profiles obtained 
in the E and NW rims with those derived from a simple Sedov solution, to derive physical 
parameters as well as to check whether or not they are well explained by a simple Sedov 
solution. Varying the value of the ambient density and calculating x^-values, we search for the 
most suitable (best-fit) ambient densities to represent the EM proffies for the E rim and the NW 
rim, respectively. Figure 5 shows the EM proffies obtained in the E (data points with circles) 
and the NW (data points with triangles) rims as well as the best-ffi EM proffies (solid line for 
the E rim and dashed line for the NW rim, respectively) as a function of radius. The best-ffi 
pre-shock ambient densities for the E and NW rims are derived to be 0.091 and 0.084 cm~^, 
respectively. These values are consistent with previous results based on X-ray observations 
(0.01-0.2 cm~^). Looking at figure 5, we see apparent deviations at the outermost regions for 
both the E and NW rims; the data exceed the best-ffi models by factors of ~1.8 (E) or ~1.4 
(NW). This fact suggests that the forward shock is not propagating into a uniform ambient 
density, but recently encountered a denser material than the shock had propagated before. The 
gradual increase of r toward the shock front (see. Table 1 or Table 2) also supports the idea 
that the forward shock is now encountering a denser material. 

Even if the SNR is not explained by a simple Sedov solution, we can estimate the age 
of the SNR by applying a simple Sedov solution in a uniform ambient density. Then, the age 
is calculated to be 11000 (Ts/0.3keV)-°-^(i?/50')(rf/lkpc) yr. Considering that the forward 
shock is now encountering a denser material, the age from the simple Sedov solution would be 
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Table 3. Spectral-fit parameters for the central ellipse and the surrounding region in the central FOV. 



Component 


Parameter 


-Cjiupse 


Surrounding 




Model: zViNEls- 


f iNEI 




Absorption 


AinixlO cm J . . . 


4.3±0.1 


4.6±0.1 




Kic[K;eVJ 


(J.o (tixedj 


0.3 (hxedj 




log(r/cm ^sec) .... 


1U.4d±U.(J5 






EM' [x 10 cm °J . . 


1 A A 1 A A r 

1.09±0.05 


A A r 1 A A O 

0.95±0.03 




Surface brightness^ 


1.3x10 


1.1x10 


VNJil (ejecta) 


kJo[kcVJ 


A AC) 1 A A 1 

0.48±0.01 


A C /i 1 A A 1 

0.54±0.01 




0(=N) 


A /I O 1 A A /I 

0.48±0.04 


2.8±0.3 




INe 


A O IT 1 A AO 

0.25±0.03 


A A 1 A 1 

0.9±0.1 




Mg 


A /I A 1 A A O 

0.40±0.03 


1.5±0.1 




Si 


1.8±0.1 


9.5±0.4 




s 


3.1±0.6 


10±1 




Fe 


r\ O O 1 AAA 

0.d6±0.02 


2.6±0.07 




log(T/cm sec) .... 


1 1 AO 1 A AO 

11.08±0.03 


11 O A 1 /'I A A 




EM' [X lU cm J . . 


A 'lO 1 A AO 

0.7d±0.02 


A 1 /I A 1 A AAO 

0.149±0.003 




Surface brightness^ 


1.2x10 


0.74x10 




Kie [keVJ 


>7 


>6 




Metal abundance . 


0.5 (hxed) 


0.5 (faxed) 




log(T/cm~'^sec) .... 


>12 


>12 




EMl'[xlOi8 cm-5] . 


0.021±0.004 


0.014±0.003 




Surface brightness^ 


2.7x10^15 


1.8x10-15 


xVd.oi 




1281/1100 


1596/1331 






Model: 2VNEIs+power-law 


Absorption 


7VH[xl02icm-2] ... 


4.2±0.1 


4.2±0.1 


VNEI (ISM) 


fcTc[kcV] 


0.3 (fixed) 


0.3 (fixed) 




log(r/cm~'^scc) .... 


10 fS2+°-°<' 

iU.OZ_o Q5 


10.651^;°^ 




EM'f[xl0iScm-5] . . 


0.93±0.04 


0.71±0.03 




Surface brightness^ 


1.1x10"" 


0.9x10-13 


VNEI (ejecta) 


fcTc[keV] 


0.50±0.01 


0.55±0.01 




0(=N) 


0.52±0.04 


2.6±0.2 




Ne 


0.27±0.03 


0.9±0.1 




Mg 


0.43±0.04 


1.5±0.1 




Si 


1.8±0.1 


8.8±0.3 




S 


2.9±0.5 


9.3±1.2 




Fe 


0.68±0.02 


2.36±0.07 




log(r/cm~'^scc) .... 


11.12±0.03 


11.32±0.02 




EMt[xlOiScm-S}5. . 


0.65±0.01 


0.15±0.003 




Surface brightness^ 


1.0x10^13 


0.69x10-13 


Power-law 


Photon Index 


1.4±0.1 


1.5±0.2 




Norm-'- 


3.9±0.7 


2.9±0.6 



3x10" 



10'* 




I ' ' I 

3x10" 3.5x10" 4x10" 

Radius (cm) 

Fig. 5. EM profiles as a function of radius. Our data are shown with circles for E rim or triangles for 
NW rim. The solid or dashed lines are the best-representative models for E rim and NW rim, respectively. 
We set the radius of the SNR to be i? = 48' (4.3x10^^ cm at a distance of 1 kpc) which is derived from 
the ROSAT PSPC surface brightness contour. 

an upper limit of the age of the SNR. 

We found possible evidence that the forward shock is now interacting with a denser 
material which surrounds the SNR. A candidate for the denser material is indeed found around 
G156.2-I-5.7 as an HI shell surrounding this SNR (Reich et al. 1992). Therefore, the pre-shock 
ambient density of ~0.09cm~^ derived by the comparison with a simple Sedov profile must be 
an upper limit for the density inside the denser material. In this context, the SN explosion which 
produced G156.2-I-5.7 occurred in a very low-density (less than 0.09 cm~^) medium. The low 
density is reasonably expected to be a cavity evacuated by a stellar wind of the progenitor star. 
On the other hand, the relatively dense material that the forward shock is now encountering is 
suggestive to a wall of the cavity. 

A powerful method to estimate the mass of the progenitor star is to compare relative 
metal abundances measured for the ejecta with those expected in nucleosynthetic models. We 
successfully measured the metal abundances of the ejecta, eliminating the contamination from 
the swept-up ISM along the line of sight by using a two- (i.e., the ISM and the ejecta) component 
VNEI model. Figure 6 shows the comparison of several metal abundances relative to O between 
our data and various nucleosynthetic models. Here, we take the metal abundances measured in 
the surrounding region in the central FOV as the representative ejecta abundances. We employ 
Type-la models from Iwamoto et al. (1999), while core-collapse models with progenitor masses 
of 11 Mq from Woosley & Weaver (1995), 13 M© from Thielemann et al. (1996), 15 Mq, 21 M©, 
and 25 M0 from Rauscher et al. (2002). We then find that our data prefer core-collapse models 
with relatively low-mass progenitors (below 15 M©) for the origin of G156.2-I-5.7. However, we 
should note that the metal abundances of the ejecta are determined in a very limited region 
of 0156.2+5.7. In order to determine the metal abundances for total ejecta, we need further 
X-ray observations of 0156.2+5.7. 
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Ne Mg Si S Fe 



Fig. 6. Metal abundances of Ne, Mg, Si, S, and Fe to O. Our data are those derived for tiie ejecta 
component at tiic surrounding region in the central FOV. Red solid lines are those expected in Type-la 
(CDDl or W7) models (Iwamoto et al. 1999). Blue lines are those expected in core-collapse models from 
relatively low-mass progenitors (11 M© [Woosley & Weaver 1995]; ISM© [Thielemann et al. 1996]; 15 M© 
[Rauschcr ct al. 2002]). Green lines also represent core-collapse models but from relatively high-mass 
progenitors (Rauscher et al. 2002). 

The HI shell surrounding G156.2-f5.7 is possible evidence of the wall of the cavity as 
already suggested by Reich et al. (1992). If this HI shell is indeed the wall of the cavity, the 
radius of the cavity is estimated to be ~15 ((i/lkpc)pc (i.e., the radius of the SNR). McKee 
et al. (1984) calculated the size of a cavity created by a massive star of type earlier than BO 
(~15Mq) to be a radius of 50 pc. Within the range of the distance (1-3 kpc) to G156.2-I-5.7, 
the radius of the cavity is lower than 50 pc. Therefore, the relatively small size of the possible 
cavity also supports that the progenitor star of G156.2-f5.7 is less massive than 15 M©. 

4-2. Non- Thermal Emission 

In the NW rim, the non-thermal emission comes from just behind the forward shock 
(see. Table 2). This suggests that the non-thermal emission originates from the relativistic 
electrons accelerated by the forward shock. Then, we should note that the value of kT^ just 
behind the forward shock is derived to be about 0.3 keV, which results in a shock velocity 
of ~500 (A;Tc/0.3 keV)°'^ kmsec"^ assuming thermal equilibration between ions and electrons. 
The assumption seems to be reasonable, since the the expected proton temperature, ~0.2 
(r/6x 10^° cm~^sec)~^(/i;Tc/0.3keV)^'^ keV, is comparable to the measured electron tempera- 
ture. Although this equation is valid only for highly non-equilibrium stage of shocked plasma, 
the resultant value at least supports that thermal equilibrium is almost reached. In addition, 
we should note that the electron-proton thermal equilibration is expected in slower shocks (e.g., 
Rakowski 2005). It is surprising that we detect non-thermal X-ray emission from such a slow 
shock, since a currently accepted theory of diffusive shock acceleration requires high shock ve- 
locity of at least 2000kmsec"^ to produce non-thermal emission in the X-ray energy band. The 
slow shock might require unacceptably high efficiency of particle acceleration. Therefore, we 
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Table 4. Spectral-fit parameters derived by the SRCUT model. 



Parameter 


NWl 


NW2 


NWS 


Ellipse 


Surrounding 


J^rolloff[xlOl6Hz] ... 


0.95±0.05 


1.28±0.04 


1.28±0.02 


1.63±0.09 


1.20±0.05 


Surface brightness^ 


0.75x10-1'' 


1.2x10-" 


1.2x10-" 


0.69x10-" 


0.43x10-" 


xVd.o.f 


691/590 


1043/910 


972/797 


1315/1101 


1610/1332 



*The errors are in the range Ax^ < 2.7 on one parameter. 
^In unit of ergcm~^sec~^ arcmin"'^ in a range of 0.5-10 keV. 



check the vahdity of non-thermal origin for the hard-tail emission, based on morphologically 
and spectral features. 

Regarding the morphology, relatively strong non-thermal radio emission is detected in 
the northern and the southern shell regions, while relatively weak emission is detected from 
the central portion of the remnant. Therefore, the distribution of the X-ray hard-tail emission 
matches the radio morphology; the X-ray hard-tail emission is detected with strong intensity 
in the NW rim, with weak intensity in the central portion, while it is not significantly detected 
in the E rim. This suggests that the radio emission and X-ray hard-tail emission have the 
same origin with each other, supporting the origin of the X-ray hard-tail emission is non- 
thermal synchrotron. In this context, we expect to detect non-thermal X-ray emission from the 
southern rim of the SNR, which will be checked by future X-ray observations. 

To check spectral features in multi-wavelength emission, we introduce the SRCUT model, 
which describes synchrotron radiation from a power-law distribution of electrons with an expo- 
nential cut-off (Reynolds & Keohane 1999) instead of the power-law model. In this model, we 
fix the radio spectral index to 0.48 (Xu et al. 2007). The flux at 1 GHz is also fixed at the value 
determined by radio observations (based on figure 1 in Xu et al. 2007); 2.1 mJy arcmin"^ for 
the NW rim, 0.83mJyarcmin~^ for the central portion. Thus, the only free parameter in the 
SRCUT model is the roll off frequency, I'mUoS- In this way, we fit all the spectra in which we 
detected hard-tail emission. The measured roll off frequencies and x^-values are summarized in 
Table 4. The fits are fairly good, although they are slightly worse than those derived by using 
the power-law model (see. Table 2 or Table 3). The measured roll off frequency, ~1.5xl0^^ Hz, 
is relatively low compared with those derived in young historical SNRs; 26xlO^^Hz for Tycho 
(Tamagawa et al. 2008), 20xl0i^Hz for Kepler (Cassam-Cheinai et al. 2004), and 5.7x10^6Hz 
for SN1006 (Bamba et al. 2008). The low value of t'roiioff is in a sense that the efficiency of 
particle acceleration can be reduced. Now, we calculate the efficiency of particle acceleration. 
The roll off frequency is expressed as 

i^rouoff = 1.6 X 10^^(5/10/iG)(E^ax/10TeV)2 Hz, 

where B is a. strength of ambient magnetic field and -Emax is a maximum energy of accelerated 
electrons (Reynolds 1998). Equating the synchrotron loss time to acceleration time, -Emax is 
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written down as follows (Uchiyama et al. 2003), 

^max = 67(5/10/iG)"°-^(i;s/2000kmsec~^)r/-^/2 TeV, 

where Vs is the shock velocity and r] is the so-called gyrofactor, i.e., the ratio of the particle mean 
free path to the gyroradius, and represents the efficiency of particle acceleration. Combining 
the above two equations, we can derive the following equation, 

i^rouoff = 7.2 X 10^^(t;s/2000kmsec"^)^r/-^ Hz. 

Substituting the measured t'roiiofr of 1.5x 10^^ Hz and the forward shock velocity of 500kmsec~^, 
we obtain the value of t] to be 48 (If thermal equilibration has not achieved yet behind the shock 
front, the shock velocity will be faster than 500kmsec~^ so that t] will be larger than 48). Thus, 
the particle acceleration is much less efficient than that taken place in the RXJ1713. 7-3946 SNR, 
in which we see the extreme Bohm diffusion limit, 77 = 1 (Uchiyama et al. 2007). Therefore, the 
fact that the slow (~500kmsec~^) shock in G156.2+5.7 produces non-thermal X-ray emission 
turns out to be acceptable from a theoretical point of view of diffusive shock acceleration. We 
consider that the detection of non-thermal X-ray emission associated with the slow shock in 
G156.2+5.7 is just simply because the flux ratio of non-thermal to thermal emission in this 
SNR is higher than those in the general middle-aged SNRs where non-thermal X-ray emission 
is not detected. Ellison et al. (2000) argued that a weak ambient magnetic fleld and/or a low 
ambient density yield the high flux ratio between non-thermal and thermal emission. As for 
the ambient density, we obtain similar ambient densities between the E rim (0.093 cm~'^) and 
the NW rim (0.087 cm~^). However, we did not detect signiflcant non-thermal emission from 
the E rim, but detected it only from the NW rim. Therefore, we suggest that the magnetic 
fleld plays an important role to the non-thermal/thermal flux ratio for the case of G156.2-I-5.7. 

We discuss the nature of the central non-thermal emission in G156.2+5.7. The central 
bump isolated from both the northern and the southern shell in the radio image reminds us 
of a plerion for the nature of the central non-thermal emission, although there is no apparent 
counterpart of a pulsar/pulsar wind nebula. Xu et al. (2007) conclude that the emission from 
the central portion came from the shell of the SNR rather than a plerion because radio spectra 
do not show spectral flattening in the central portion of the remnant which is expected in the 
case of pulsar wind nebula. In principle, it is possible for the shell to be projected to the 
center of the SNR. In the simulations of Orlando et al. (2007) to investigate the origin of the 
asymmetries in bilateral SNRs, we flnd such an example. In the paper, we flnd that the radio 
morphology of G156.2+5.7 is similar to the surface brightness of non-thermal emission expected 
in a situation that the bilateral radio limbs are parallel to the ambient magnetic fleld and a 
gradient of ambient density is perpendicular to the magnetic fleld (see flgure 5 in Orlando et 
al. 2007). Although the possibility that the central non-thermal emission is a plerion can not 
be completely ruled out, it is preferred that the central non-thermal emission is associated with 
the shell of the SNR. Further X-ray observations around the central portion of the SNR will 
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allow us to investigate a detailed variation of spectral index around the center of the SNR. If 
we see strong spectral variation, as is generally seen in pulsar wind nebulae (e.g., Petre et al. 
2007), the origin would be most likely a plerion, otherwise the central non-thermal emission is 
associated with the shell of the SNR. 

5. Conclusions 

We have observed the middle-aged SNR G156.2+5.7 with the Suzaku satellite in three 
pointings which are targeted at the E rim and the NW rim, and the central portion of the 
SNR. Hard-tail X-ray emission is detected in the NW rim and the central portion, while it is 
not significantly detected in the E rim. The distribution of the hard-tail emission seems to 
be coincident with that of the radio non-thermal emission. Also, multi-wavelength emission is 
well represented by a so-called SRCUT model. These facts lead us to conclude that the hard- 
tail emission originates from non-thermal synchrotron emission. The velocity of the forward 
shock which likely produces the non-thermal X-ray emission in the NW rim is estimated to be 
~500kmsec~^, assuming thermal equilibrium between ions and electrons. The shock might be 
the slowest among those having non-thermal emission in the X-ray wavelength. The relative 
abundances in the ejecta component suggests that the origin of G156.2+5.7 is a core-collapse 
SN explosion whose progenitor is a relatively low-mass star of ~15M0. 
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